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Both three- and one-dimensional studies of local multiplicity fluctuations in 
hadronic Z decay are performed using data of the L3 experiment at LEP. The 
normalized factorial moments in three dimensions exhibit power-law scaling, indi¬ 
cating that the fluctuations are isotropic, which correspends to a self-similar fractal. 

A detailed study of the corresponding one-dimensional moments confirms this con¬ 
clusion. However, two-jet subsamples have anisotropic fluctuations, correspending 
to a self-affine fractal. These features are, at least qualitatively, reproduced by the 
Monte Carlo models jetset and HERWIG. 


1 Introduction 


Dynamical fluctuations can be investigated using the normalized factorial mo¬ 
ments (NFM), 
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where M is the number of bins in which momentum space is partitioned and 
n m is the multiplicity in the mth bin. If power-law scaling (intermittency), 

F q (M) <x M+' , (2) 


is observed then the corresponding hadronic system is fractal, as is ex¬ 
pected from a branching process. Since particle production occurs in three- 
dimensional momentum space, we can distinguish two cases: If the scaling is 
observed when the space is partitioned equally (unequally) in each direction 
then the fractal is self-similar (self-affine). For two dimensions the degree of 
anisotropy is quantified by the so-called Hurst exponent, H a b = In M a / In Mb, 
where M a and Mb are the number of bins in the two dimensions. Projected 
onto one dimension, a, the second-order NFM saturates and is given byld 

F^ a \M a ) = A a - B a M~^ a . (3) 
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H ab is related^ to the exponents 7 
through H ab = (1 + 7 b)/(l + 7a)- 
In three dimensions, the dynamical 
fluctuations are isotropic if H ab = 1 
for all pairs of axes, ab, or, equiva¬ 
lently if the 7 a for all three axes, a, 
are equal. 

Consequently, we measure the 
NFM’s in three dimensions to exam¬ 
ine whether there is power-law scal¬ 
ing and to measure the intermittency 
indices, cf> q . We also measure the 
one-dimensional F% to get the expo¬ 
nents 7 a for the three directions. 

To investigate possible depen¬ 
dence of the fluctuations on pertur¬ 
bative or non-perturbative scales, we 
also examine the F 2 for 2 -jet events 
as a function of y cut . 

2 Method of analysis 

Since we are studying dynamical 
fluctuations, the coordinate system 
should not depend on the dynamics 
of qq —> hadrons. So we choose a 
cylindrical coordinate system with z 
axis along the qq direction and we 
partition the phase space in bins of 
y 1 pt, p. The y and pt directions are 
defined by the electro-weak process. 
But, we need an origin for ip. Clearly, 
we should not use an axis, like major, 
which depends on the QCD dynam¬ 
ics. One possibility!] is to rotate the 
coordinate system around the thrust 
axis and to define a new x axis as ly¬ 
ing in the z-bearn plane. However, it 
is easier simply to choose the origin 
at random. 
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Figure 1. The 3D NFM as a function of 
the number of bins, M, compared with re¬ 
sults of JETSET and HERWIG. 
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Figure 2. The F 2 in one dimension as 
function of number of bins, M, in the ran¬ 
domized frame, compared to jetset and 
HERWIG. 
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Table 1. The fit parameters of the three-dimensional NFM 



Randomized 

Corrected to qq 

order 

09 

xV dof 

09 

xVdof 

2 

0.194 

± 

0.003 

± 

0.003 

8/9 

0.221 

± 

0.003 

± 

0.003 

8/9 

3 

0.598 

± 

0.011 

± 

0.014 

11/9 

0.685 

± 

0.011 

± 

0.012 

11/9 

4 

1.082 

± 

0.013 

± 

0.018 

8/9 

1.206 

± 

0.022 

± 

0.026 

7/9 

5 

1.731 

± 

0.024 

± 

0.025 

6/9 

1.858 

± 

0.028 

± 

0.035 

7/9 


Experimentally, we estimate the qq direction by the thrust axis. We also 
examine the effect of correcting this direction to the qq direction using a 
correction factor determined from JETSET. 

The data used in the analysis were collected by the L3 detector in 1994 
at the Z pole. Hadronic events are selected and cuts are applied to obtain 
well measured tracks and calorimetric clusters. A total of about one million 
events satisfy the selection criteria. 

The NFM’s are calculated from the raw data and corrected for detector 
effects by a factor determined from generator level MC and detector level 
MC. Systematic errors are estimated from different selection cuts, different 
methods of event selection and different MC for detector correction. 

3 Results of F q for the full sample 

The three-dimensional F q for an isotropic partitioning of phase space, are 
shown in Fig. [j] together with the results of fits of Eq. |], omitting the first 
point to eliminate the influence of momentum conservation. The results of 
the fit are also given in Table[I]. The fits are good, as is expected if dynamical 
fluctuations are isotropic. Correcting from the thrust to the qq direction 
systematically increases the values of the 4> q . Both jetset and HERWIG are 
seen to agree well with the data. 

The one-dimensional F 2 are plotted in Fig. || and the results of fits of Eq. |] 
are shown in Table|^. The values of 7 for the three axes are equal, confirming 
that the fluctuations are isotropic. Correcting to the qq direction decreases 
slightly the 7 values. Again, jetset and herwig agree with the data. 

4 Results from 2-jet subsamples 

We use the Durham algorithm to define 2-jet events for various values of the 
jet resolution parameter, y c u t- We fit the resulting one-dimensional F 2 using 
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Eq.|^ and plot the values of 7 vs. kt = sj sy cu t in Fig.[|. While 7 increases 
with k t (or y cu t) for each of the variables, the dependence on y cut is very 
different for the different variables, y, pt , and ip. JETSET and herwig behave 
qualitatively similarly to the data. 

5 Conclusions 

In three dimensions, the factorial moments exhibit power-law scaling when 
momentum space is partitioned isotropically. Fits to the three-dimensional 
F q vs. M determine the intermittency indices (j> q . In one dimension, fits to 
F 2 (y), F 2 (pt), F’li'p) confirm the isotropy. Both JETSET and herwig agree 
remarkably well with both the three-dimensional F q and the one-dimensional 
F 2 . For 2-jet events, F 2 (y ), F 2 (p t ), F 2 (p) depend very differently on the jet 
resolution parameter. Similar behavior is seen in JETSET and HERWIG. 


Table 2. The fit parameters of the 1-D 
NFM. 


Randomized 


7 

X^/dof 

y 

0.992 ±0.014 ±0.028 

36/36 

Pt 

0.986 ± 0.022 ±0.032 

29/37 

P 

0.993 ±0.024 ±0.033 

25/35 

Corrected to qq 


7 

X 2 /d°f 

V 

0.966 ±0.019 ±0.018 

35/36 

Pt 

0.972 ± 0.042 ±0.030 

31/37 

<P 

0.967 ±0.034 ±0.024 

27/35 
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